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Summary 

This  paper  contains  an  exact  solution  in  series  fora  to  the  problem 
presented  b.v  the  irrotational,  axisymraetric  flow  of  an  ideal,  incompressible 
fluid  past  a solid  torus  of  circular  cross  section.  At  infinity  the  fluid 
is  assumed  to  be  in  a state  of  unifora  motion  parallel  to  the  axis  of  the 
torus*  The  solution  is  based  on  the  use  of  toroidal  coordinates,  end  is 
given  in  terms  of  Legendre  functions  of  fractional  order  as  'em  11  as  com- 
plete and  incomplete  elliptic  integrals  of  the  first  and  second  kind*  The 
individual  co^-su«ut  solutions  employed  here  are  interpreted  physically 
through  the  j*  r^latlou  to  the  basic  ring  singularities  represented  by  the 
source  ring  and  the  vortex  ring. 

Ths  problem  is  first  approached  via  Stokes*  stream  function  exclusively 
and  is  subsequently  received  independently  in  terms  of  the  velocity  po- 
tential alons*  The  convergence  of  the  pertinent  series  developments  is 
found  to  be  unusually  favorable,  end  s complete  stream  line  pattern,  car- 
responding  to  an  illustrative  numerical  example,  is  included* 
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Introduction 


Potential  problems  involving  toroidal  boundaries  have  received  re- 
peated attention  in  the  literature.  Apparently  the  first  investigation 
of  this  kind  is  duo  to  Carl  Neumann  Qlj^  sho  dealt  with  the  distribution 
of  temperature  in  a solid  torus  and  with  the  analogous  problem  in  electro- 
statics. The  corresponding  axisymetric  problem  in  ideal  fluid  flow  was 
first  investigated  by  Hicks  £2^  and  was  reconsidered  later  by  Dyson  £3^» 
who  also  examined  certain  unsywaetric  flow  canes  associated  with  the 
torus.  The  rigorous  treatment  of  allied  problems  in  the  theory  of  elasti- 
city is  more  recent;  Freiberger  disposed  of  the  problem  of  pure  torsion 
of  a torus  of  circular  cross  section,  whereas  the  case  of  purs  bending  was 
treated  by  the  current  authors  £ff].  It  is  this  last  investigation  which 
motivated  the  present  study,  and  the  analytical  prerequisites  needed  here 
were  established  to  s considerable  extent  In  connection  with  the  work 
described  in  [jf]. 

Although  the  hydrodynamic  investigations  cited  previously  yielded 
formal  solutions  characterising  the  lmtational  axisymmtrlc  flow  of  an 
inviscid,  incompressible  fluid  past  a torus  of  circular  centerline  and 
circular  cross  section,  the  problem  appears  to  be  in  need  of  further  at- 
tention. Bloks  £2]]  based  his  analysis  on  toroidal  coordinates  and  toroidal 
harmonics  but  did  not  deal  with  the  determination  of  the  stream  function, 
which  is  of  primary  physical  interest.  Dyson  £3 j],  on  the  other  hand,  arrived 
at  a representation  of  the  potential  and  stream  functions  in  terns  of  series 
of  definite  integrals;  the  significance  of  this  representation  for  numeri- 
cal evaluations  of  the  results  appears  to  be  rather  limited.  In  view  of 

^Bashers  in  brackets  refer  to  the  bibliography  at  the  end  of  the  paper. 
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the  currently  available  nusarlc&l  tabulations  of  the  Legendre  functions  of 
fractional  order  [VQ;  a re-examination  of  the  problem  would  seem  to  serve 
a useful  purpose* 

In  19U8  Weinstein  [_*7[]  discussad  axially  syasetric  potential  flow3  and, 
in  particular,  considered  an  extension  of  the  me'  -d  of  sources  and  sinks 
to  cases  in  which  the  singularities  are  no  longer  confined  to  the  axis  of 
symmetry.  This  paper  stimulated  further  interest  in  related  problems.  Thus, 
Streeter  QT],  with  the  aid  of  ring  and  disk  singularities,  obtained  stream- 
line patterns  corresponding  to  flows  around  various  blunt-nosed  half-bodies 
and  annular- aha pod  bodies;  his  approach  is  based  on  a seal-numerical  pro- 
cedure. At  the  same  time  elliptic  integral  representations  were  established 
[V]  for  the  potential  and  stream  functions  appropriate  to  rotationally 
nystric  distributions  of  sources  and  vorticity  over  circular  rings  and 
disks.  Concurrently,  Van  Tuyl  deduced  a similar  analytical  repre- 
sentation of  the  flow  generated  by  a uniform  source  disk,  and  applied  it 
to  the  study  of  a new  family  of  half-bodias.  Still  more  recently, 

Streeter  arrived  at  the  flow  t round  a torus  of  nearly  circular  cross 
section  by  superposition  vf  a ring  doublet  upon  a uniform  stream.  This 
solution  will  later  be  identified  as  a first  approximation  to  the  rigorous 
solution  for  the  torus  of  exactly  circular  cross  section;  moreover,  tho 
results  to  be  developed  here  are  free  from  the  indirect  trial-and-error 
features  inherent  in  Streeter's  approach.  Reference  should  also  be  made 
to  a paper  by  Schiffmau  and  Spencer  £l2j  in  which  toroidal  coordinates  are 
used  to  study  the  flow  aooui  - lens-shaped  object* 

In  what  follows  we  examine  certain  pertinent  sequences  cf  singular  so- 
lutions of  the  governing  equations  referred  to  torcide.1  coordinates.  The 
physical  significance  of  these  component  solutions,  in  terms  of  *hich  the 


solution  to  the  torus  problem  will  ultimately  be  expanded,  is  investigated 
with  particular  attention  to  the  presence  or  absence  of  a circulation  in 
the  corresponding  potential  and  stream  functions:  it  is  shown  how  the  fore- 

going component  solutions  may  be  generated  through  s\:ccessive  limit  pro- 
cesses applied  to  the  basic  axisymmetric  ring  singularities  of  the  source 
ring  and  the  vortex  ring.  This  part  of  the  paper  may  be  of  interest  beyond 
the  present  application.  The  torus  problem  is  then  formulated  and  solved 
on  the  basis  of  the  stream-function  approach.  A second,  independent  so- 
lution, roa ting  on  the  potential-function  approach,  1*  indicated  merely 
since  it  provides  further  insight  into  the  problem  under  consideration. 

In  order  to  render  the  present  exposition  sensibly  self-contained,  some  of 
the  material  given  in  £5^  has  also  been  briefly  included  here. 

The  governing  Equations.  Toroidal  Coordinates 

For  the  sake  of  convenience  we  suarw’-iie  at  this  place  the  basic 
equations  governing  the  steady  irrotational  flow  of  an  ideal  incompressible 
fluid  in  the  presence  of  axial  symmetry.  Vlth  reference  to  the  cylindri- 
cal coordinates  (p3y,z),  where  the  z-axis  is  assumed  coincident  with  the 
axis  of  symnetry,  the  velocity  potential  0(p,s)  satisfies  Laplace's 

* /Mi  • ♦ 4 ae  A n 

w ^ w<—  v — t...  ■ « ■ — v 4 

* mW  ■ 0 <l>2 

and  the  velocity  components  are  given  by 

V = 0 , v = 0 (2) 

p *p*  z v ' 

Stoke:-'  stream  function  $^{p,z)  is  introduced  through  the  relations, 
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Subscripts  attached  to  functions  which  originally  bear  no  subscript 
denote  differentiation  with  respect  to  the  argument  indicated* 
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(3) 


and,  alternatively,  admit3  the  representation, 

^(f,i)=  | p0z  dp  - p0p  ds J 

<Po'‘o> 


CU> 


in  which  (p0,zQ)  are  the  coordinates  of  an  arbitrary  fixed  point  In  the 
meridional  half -plane  p - 0.  In  view  of  Equation  (1),  which  may  be  re- 
garded as  an  integrability  condition  for  Equations  (3),  the  line  integral 
(U)  is  independent  of  the  path,  and  ^(p,s)  is  single-valued,  in  any  simply 
connected  domain  throughout  which  0(p,s)  is  regular.  As  was  emphasized 
by  Veinstein  the  single-valuedness  of  V can  no  longer  be  inferred 

if  the  domain  of  regularity  of  0 is  multiply  connected.  Similarly,  "xud 
under  analogous  circumstances,  a single-valued  stream  function  nay  give 
rise  to  a multiple-valued  velocity  potential.  A potential  function  and  a 
stream  function  which  are  related  according  to  Equations  (3)  will  hence- 
forth be  called  conjugate  with  respect  to  each  other. 

It  follows  from  Equations  (1),  (3)  that  conforms  to  the  "stream 
equation". 


h h fc)  = 0 

OO  O p u2  p l 


and  Equations  (2),  (3)  further  imply. 


(5) 


(6) 


A specific  velocity  field  determines  its  associated  potential  as  well  as 
the  conjugate  stream  function  only  within  an  arbitrary  additive  constant. 
Finally,  we  recall  that  the  curves  ^(p,»)  = constant  are  the  stream 
lines  and  that  2 ~ constitutes  the  total  flow  between  the 

stream  surfaces  generated  by  ifj-  and  tyJ- 


In  what  follows  we  chill  have  occasion  to  refer  to  toroidal  coordinates 


(a,p,  y)  which  are  defined  through  the  transformation, 


where  (x,y,z)  are  the  Cartesian  coordinates  and 

p = cos  c<»  p = sin  a,  q - cosh  p,  q = sinh  p (8) 

The  ranges  of  the  curvilinear  coordinates  are  given  by 

o - d-  2rr , o * p < oo,  o - y-  2rr  (9) 

so  that 

- 1 - p,  p - 1,  i - q < go,  0 - q < co  (10) 

Furthermore,  with  the  auxiliary  notation, 

1/2 

P = (q  - P)  (11) 

we  have, 

P » (*2  ♦ y2)V2  « -5*  02) 

^2 

Equations  (7),  (8),  (11)  imply 

(p  - §)  + ^ p^  + (»  **  €)  = jrj  (13) 

The  coordinate  surfaces  x=  constant  are  thus  spherical  bowls  of  radius 

| p | centered  on  the  z-axlc  at  z 17  p/p.  The  surfaces  p - constant 

are  tori,  the  meridional  cross  sections  of  which  are  the  circles  of 
_-l 

radius  q centered  on  the  p-axia  at  p = q/q.  As  p — ► co  the  tori 
degenerate  into  the  circle^  p = 1,  i = 0,  Figure  1 shows  the  traces  of 

3 -----  - 

This  normalization  of  the  coordinate  transformation  (7)  is,  of  course, 
unessential. 
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the  foregoing  two  rutually  orthogonal  families  of  coordinate  surfaces  upon 
the  meridional  half-planes  y — constant. 

The  differential  of  arc-length  is  defined  by 

d.*  • ♦ (|&)2  ♦ (&  (Ik) 

h h2  h3 

and  the  metric  coefficients  here  appear  as 

9 2 

\ «=  h2  = fj}t  h3  = A-  (15) 

q 

A routine  computation  now  yields  for  Laplace's  Equation  (1)  in  toroidal 
coordinates, 

“ P^a+  (1  ‘ = 0 (16) 

whereas  the  stream  Equation  (5)  becomes 


^ ^p)  ♦ pq  (1  - pq)  ^p=  0 (17) 


The  toroidal  transform  of  Equations  (3)  is 


*ri  k 


(18) 


and  the  toroidal  components  of  the  velocity  vector  assume  the  form. 


h-x=- 4 % 

q “ 


Y (19) 


J 


Toroidal  Potential  and  Stream  Functions^ 

In  this  section  we  consider  certain  aggregates  of  solutions  to  La- 
place's Bquation  (16)  and  to  the  stream  equation  (17),  in  toroidal  coordi- 

See  £l33>  Chapter  X,  for  a discussion  of  toroidal  harmonics.  Toroidal 
stream  functions  are  also  dealt  with  in  [”5^« 
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nates.  Neither  of  these  tiro  equations  Is  separable  In  the  strict  sense. 
There  do  exist,  however,  pseudo-product  solutions  of  the  form, 


0(ot,p)  * yU.A(a)B(p) 


(20) 


to  Laplace's  equation,  whereas  the  stream  equation  admits  the  pseudo-product 
solutions. 


if 


C(a)D(^) 


(21) 


Substitution  of  Equations  (20),  (21)  into  Equations  (16),  (17),  subsequent 
separation  of  variables,  and  adherence  to  the  restriction  that  the  solutions 
sought  have  a period  2 TT  in  cl,  ultimately  yield, 


~ /*£cos  na  or  sin  naj  £pn_  or  Si-l/2^J 

^(»^)  = jj^coa  not’or  oUi  *“*]  [PA-l/2^  °r  Qn-l/2(q^J 

(n  — 0,  1,  2,  • . • ) 


(22)- 


where  Qn  ^ tr®  the  I*8®ndh'®  functions  of  the  first  and 

second  kind,  respectively. 

For  future  reference,  we  cite  Legendre ' s equation,  ■ - 


+ - (-2  - E)pn-  1/2  * ° 


(23) 


and  recall  the  recursion  formulas. 


-n  - 1/2  Pn-l/2 


**“*0-1/2  “ (2n+  1)pn+  l/2  + (2n’1)Pn-3/2 
2q2pn-l/2  " (2n”D(qPn-l/2  " Pn-3/2^ 


(2U) 


?The  primes  denote  differentiation  with  respect  to  the  argument  q. 


Equations  (2U)  are  valid  for  all  integral  values  of  n,  and  remain  valid 


if  Pn  jyg  is  replaced  with  Qn  Tie  also  note  the  bilinear  identity, 


Pr.  - l/2*h  - 3/2  ~ PA  - 3/2Qn  - 1/2  = ” 

Furthermore,  we  shall  need  to  refer  to  Laplace's  integral, 

r * 


P-l/2(q)  " W 


da 


+ a 


q + q cos  ct 


and  to  the  Integral  representation. 


9n-  " y? 


7T 


cos  tut 


do. 


which  gives  rise  to  the  Fourier  expansion. 


(25)' 


(26)' 


(27  )* 


00 


z<  ~r  ^-i/2 T ir  s - 9n- i/j  c°*  ” 


(26)- 


where  is  defined  in  Equation  (11) • 

Equations  (26),  (27)  supply  the  link  to  the  representation  of  the  Le- 
gendre functions  of  fractional  index  in  terms  of  elliptic  integrals.  Thus, 
let 


=7^ 


cosh  p/£* 


- “ rH 


(29) 


Setting  6 = a/2  in  Equation  (26),  we  arrive  at 

rt/2 


-1/2 


2k  d®_ 

^ 4 Vl  - k2 


2 7~ 

sin  0 


2k 

■71' 


(30) 


°See  [13],  art.  U5. 

7See  [13],  p.  381. 

®See  [14],  Pc  U43* 

o 

The  argune,*t  of  Pn_i/2  ^n-1/2  henceforth  assumed  to  be 

q,  unless  otherwise  specified. 


and  setting  9 “ {it  - a)/2  in  Equation  (27),  we  reach, 


(31) 


where  K and  K'  are  the  coaplete  elliptic  Integrals  of  the  first  kind 
'referred  to  the  modulus  k and  the  conplementary  modulus  k',  respectively. 
Alternatively,  we  may  define  moduli. 


(32) 


(33) 


The  near  moduli  are  related  to  k and  k' , introduced  in  Equation  (29), 
through  the  Landen  transformation,*® 


“L 


1 - k* 

mF 


(3U) 


so  that 


TL  * (1  k' )K'  « K* 

VEJ 

B (k*  - 1)E*  + *» 


j 


Here  and  E^  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind  based  on  k^.  The  moduli  k^,  kj^  are  identical  with  those 

*°See  £li*^,  p*  507 • 

**The  second  of  Equations  (35)  is  obtained  by  differentiating  the 
first  with  respect  to  k». 


previously  employed  by  Hicks  £*j,  Ecuquet  £li>3»  the  current  authors 
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Symmetric  and  Antisymmetric  Flow  Fields 

Among  the  potential  and  stream  functions  Hated  in  Equations  (22), 
those  involving  the  Legendre  functions  of  the  second  kind  become  infinite 
as  q -w  1,  and  therefore  violate  the  regularity  requirements  on  the  z-axio 
inherent  in  the  torus  problem  to  be  considered  subsequently.  He  shall  thus 
limit  the  present  discussion  to  the  potantisl  and  stream  functions  involv- 
ing ?a_  j/2  aDC*  its  derivative,  which  remain  analytic  for  every  finite  q 
and  vanish  at  Carte elan  infinity!  they  necou>  singular  on  the  "limit  circle", 
as  q — ► oo. 

We  may  distinguish  between  two  basic  flow  cases  according  as  the 
velocity  field  is  symmetric  or  antisymmetric  with  respect  to  the  plane 
a = 0,  that  is,  according  as  the  velocity  potential  is  an  even  or  odd 
function  of  s and  a.  Symmetric  flows  are  generated  b-v  the  potentials 
and  stream  functions, 


n C MPn-l/2  008  0(1 


-2 

TJ  = pt  . sin  not 
* n fi  n-  1/2 


■ (36) 


(n  = 0,  1,  2,  ...)  J 

It  is  important  to  observe  that  and  Vn  do  not  constitute  a conju- 

gate pair  in  the  sense  of  Equations  (18).  We  now  seek  to  establish  the 
stream  function  Vn  and  the  velocity  potential  which  are  conjugate 

to  $n  and  Vn,  respectively.  To  this  end  we  confirm  with  the  aid  of 


Bquations  (16)  that, 


12 


Vn*l-^  = ^<Vn  + 1-  Vn) 


a . a 

*n  + l *n 


2n  + 1 


<*»♦!-  *»> 


(n  R 0,  1,  2,  « » « ) 


Since  VQ  = 0,  ee  may  take  $0  to  be  zero,  end  the  foregoing  recursion 


relations  yield 


vn - v «► 2 

n o 


*r%v  <wi  - *i-i> 


$n'j>  (2i  - Ut&i  - 4>t.  x) 


(n  — 1,  2,  3,  •••) 


In  order  to  conplete  the  detersinatlon  of  for  n = 0,  1,  2,  ...» 

ee  need  to  find  For  this  purpose  ue  introduce  the  pa  rare  ter  t 

through 


»(t,k)  = 0 


7T  - a 


or,  equivalently,  through 


t * F(0,k)  = 


(1  - c sin*  *)J 


am(t,k)  being  the  anplitude  function  and  F(0,k)  the  in coup  lets  elliptic 
integral  of  the  first  kind,  both  referred  to  t je  modulus  k given  in 
Equations  (29).  Thus, 


snt  * sin  0 = cos  {?,  cnt  - cos  0 = sin  »,  dnt  = (hi) 

? z y? 


i 


} 


13 


12 

where  on,  cn,  dn  denote  the  Jacobian  elliptic  functions  for  the  modu- 
lus k. 

According  to  Equations  (30),  (36),  and  (41),  we  may  write 


= H-  f-i/2  = ^ v 


Equations  (18),  applied  to  the  conjugate  pair 
Equation  (42),  now  leaa  to 

£-*[*-*■  jg 


(1*2  > 

, by  virtue  of 


aW 


o _ yg  k*k 

7 r 


' snt  cnt 

K r~ 

<hrt 


(43) 


where  E'  is  the  cosplete  elliptic  integral  of  the  second  kind  based  on 
the  modulus  k'.  Equations  (43)  are  explicitly  integrable  by  recourse  to 
the  identities  listed  on  p.  5l6  of  [l4].  There  results 


y0  s - E')t  + k2K'  - K'E(t)J 


(W») 


13 


In  which  E(t)  designates  the  incomplete  elliptic  integral  of  the  second 
kind  for  the  modulus  k*  We  note  from  Equations  (38)  that  W n and 
which  are  conjugate  to  the  pseudo-product  solutions  $>n  and  of  the 

potential  and  stream  equations,  no  longer  possess  the  pseudo-product  struc- 
ture. Indeed,  a single  potential  function  gives  rise  to  a conjugate 

stream  function  ^ ^ whose  representation  requires  an  infinite  series  of 
stream  functions  in  pseudo-product  form;  the  closed  representation  given  by 


ee  [14] 
grala. 


, Chapter  XXII,  for  a discussion  of  Jacobian  elliptic  functions 


^See 
and  lnte 

13 

It  should  be  noted  that  the  Jacobian  elliptic  functions  appearing 
here  may  be  eliminated  by  means  of  Equations  (41)  • 


) 


111 


Equations  (38),  (Ui)  was  made  possible  only  through  the  introduction  of 
elliptic  functions  and  integrals. 

Since 


an(t  + 2K)  = ant  t-  jt 
E(t  * 2X)  *=  E(t)  + 2E 
and  recalling  Legendre's  identity, 

EE'  + E'S  - EX'  = y* 


(U5) 


(U6) 


as  well  as  the  periodicity  properties  of  the  Jacobian  elliptic  functions, 
we  conclude  from  Equations  (UI4),  (38)  that 


VB(*  + 27 r,p)  - Vn(a,p)  * 2 V? 


(n  = 0,  1,  2,  ...) 


(U7) 


The  single-valued  potentials  $n  thus  correspond  to  associated  conjugate 
stream  functions  VQ  which  are  many-valued  for  - 00  < cl  < 00,  or  else 
discontinuous  along  the  cut  a**  0 in  the  meridional  half -plane  correspond- 
ing to  0 - ct-  27f  f 0 - p < 00.  This  result  is  consistent  with  our 
preriouc  observations  regarding  the  singular  character  of  as  q — *■  00. 

Or.  the  other  hand,  the  conjugate  sequences  Vn|  (n  = 1,  2,  3*  •••) 

are  both  single-valued* 

The  treatment  of  the  ant i-syme trie  flow  fields  generated  by 
*n~  I1  ,ln  na 

K “ M 1/2  c”  “ 


0.8) 


ia  strictly  analogous  to  the  preceding  discussion  of  symmetric  fields,  and 
we  may  confine  ourselves  to  stating  the  pertinent  results. 


and 


0 , (g/ 

*n*  n 


are  conjugate  pairs  of  potential  and  stream  functions,  then 
_ -2 


jrn  ( K* i * K? 

in+l-in 


(n  — 0,  1,  2,  • • • ) 


and  consequently. 


*'.--*XsrVr«ri  - 


. W-1K0, 


(n  s 1,  2,  3,  • ••  ) 


(U9) 


(50) 


Here  we  may  take  W Q~  0 since  0Q  vanishes,  whereas  0q  turns  out  to 
be  given  by 


0O  * [(*'  - K')t  - K'B(t)J  (51) 

The  functions  0^,  and  are  periodic  with  period  27 T in  at*  On 

the  other  hand,  we  find  that 


0n(tt+  27f,0)  * 0n(<t^3)  = V? 

(n  = 0,  1,  2,  ...) 

which  indicates  a circulation  in  the  velocity  potentials  0Q 
the  single-valued  stream  functions  (^n* 


(52) 


conjugate  to 


We  now  turn  to  a physical  interpretation  of  the  potential  and  stream 
functions  discussed  in  this  section.  With  regard  to  symmetric  flows,  we 
note  on  the  basis  of  Equation  (3.9)  of  Reference  [V]  that  the  potential  of 
a uniform  source  ring  of  total  strength  m and  unit  radius  (coincident 
with  the  limit  circle  cf  the  toroidal  coordinate  system)  in  our  present 
notation  appears  as 

' - w ^ <»)■ 

with  and  defined  as  in  Equations  (33),  (35).  A trivial  compu- 
tation, involving  Equations  (53),  (35),  (33),  (32),  and  (30),  identifies 


0 « W = Wa  (51.) 

V?  V? 

as  the  potential  and  stream  -tions  appropriate  to  the  source  ring  under 
consideration.  Moreover,  we  confirm1*4  with  the  aid  of  Equation  (U7)  that 
the  stream  function  possesses  a cyclic  period  of  2m,  as  is  consistent  with 
hydrodynamic  theory. 

In  a similar  manner,  guided  by  Equation  (7.5)  of  Reference  Q93»  110 
recognise 


(55) 


as  the  potential  and  stream  functions  belonging  to  the  antisymmetric  flow 
generated  by  a uniiorm  vortex  ring  of  unit  radius  and  total  circulation  F , 
The  circulation  constant  T of  the  velocity  potential  is  readily  veri- 
fied by  means  of  Equation  ( T>? ) • 

Equations  (5U),  (55)  establish  the  significance  of  the  initial  members 
in  the  sequences  of  symmetric  and  antisymmetric  flows  discussed  previously. 


^^See  Equation  (5.7). 


The  remaining  members  in  these  two  aggregates  of  singular  flows,  which  are 
characterized  by  singularities  of  progressively  higher  order  on  the  Unit 
circle,  may  be  generated  through  successive  z-differentiations  applied  to 
the  basic  singular  solutions  which  represent  the  source  ring  and  the  vortex 
ring.  Indeed,  a somewhat  tedious  computation,  depending  on  the  repeated 
use  of  the  Legendre  relations  (2b),  yields  the  result. 


e£> 

sr 


n _ (2n  - 1)  2n  ♦ 1 m 

A-TT ‘ *n-  1 + ”^n T~  *n+  1 


3 (S'. 


n _ 2n 


(56) 
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(n  — 0,  1,  2,  ... ) 


Moreover,  these  recursion  formulas  remain  valid  if  $ and  0,  as  well  as 
V and  (£/,  are  interchanged  while  the  signs  of  the  right-hand  members  are 
reversed.  As  an  illustrative  example  we  determine  the  solution  correspond- 
ing to  a ring-doublet,  whose  potential  is  the  first  derivative  with  respect 
to  z of  the  potential  far  a source  ring.  By  Equations  (5b),  (56),  (50), 
the  ring-doublet  is  thus  represented  by 


(57) 


and  these  formulas  are  readily  transformed  into  the  expressions  derived  by 
Streeter  £lll. 


^Observe  that 

(2b),  (36),  (b8). 


- and  - - ^7, , according  to  Equations 


Solution  of  the  Torus  Problem?  The  St  rears- Function  Approach 

We  are  now  in  a position  to  consider  tha  main  problem  of  this  Investi- 
gation. Let  the  boundary  of  the  torus  coincide  with  the  parametric  sur- 
face P = P0  (&J®  Figure  i).  Then,  according  to  Equations  (13), 

A = | = q0  - co.h  p0  (S8)16 


where  a and  b are  the  radii  of  the  cross  section  and  of  the  centerline, 
respectively,  and  A will  be  called  the  "shape-ratio"  of  the  torus. 

The  problem  consists  in  determining  a stream  function  (JJ  (c c£)  which 
satisfies  Equation  (17)  throughout  the  domain  0 £ a , 0 ^ p - pQ,  and 
conforms  to  the  boundary  and  regularity  conditions, 

P'ta.pj  = c,  <('(«,(»  — »•  (a,p)  — (0,0)  (59) 


where  c is  the  Initially  unknown  value  of  (JJ  on  the  surface  of  the  torus, 
and  P ^ designates  the  stream  function  associated  with  the  uniform  velo- 
city field  Vp  » 0,  v#  = 1 at  infinity* 

Guided  by  the  antisymmetric  nature  of  the  uniform  velocity  field  at 


Infinity,  we  assume  the  desired  solution  in  the  form, 

ao 


n *=  0 


(60) 


with  (JJ^  defined  by  the  second  of  Equations  (US),  and  note  that  (JJ  (ot,p) 
automatically  conforms  to  the  condition  at  infinity  since  (//£ 0,0)  = 0* 

The  coefficients  of  superposition  aQ  (n  = 0,  1,  2,  •••)  are  to  be  deter- 
mined consistent  with  the  boundary  conditions  for  p = PQ.  By  virtue  of 
Equations  (60),  (59)>  (U8),  and  (12),  this  leads  to 


^The  subscript,  zero, 
indicate  its  evaluation  for 


attached  to  any  function  of 

P = Po- 


p,  will  henceforth 
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i •* 


CO 


> ' VA-i/s^o’ co" 

n = 0 


3*- 


2/“-o 


(61) 


i 


In  order  to  expand  the  right-hand  member  of  Equation  (61)  in  a Fourier 

series,  we  recall  the  expansion  (28),  With  the  aid  of  this  expansion,  and 

making  use  of  the  identities  (2h)  applied  to  the  Legendre  functions  of  the 

second  kind,  we  thus  arrive  at 
oo 

Vi-  V2(4o)  cu* 

coo  net 


n 


2^yr  *£_* 

~ ^ - 1 J 1 + 8no 


(62) 


in  which  8 designates  the  Kroneckar  delta.  Equation  (62)  implies 
no 


_ , 2c 

= (— j 


lin  - 1 


*n 


(63) 


provided 


knr- 


(6k) 


LK-i/M 

and,  in  accordance  with  Equation  (60),  the  appropriate  velocity  potential 
appears  as 


oo 


lf(<x*p)  = * + «n  tr 

n = 6 


(65) 


the  functions  0Q  being  defined  in  Equations  (50)  and  (51)« 


It  is  apparent  fro*  the  preceding  developments  that  the  boundary  con- 
ditions (59)  do  not  suffice  to  date  rains  the  yet  unknown  value  c of  Ip 
on  the  surface  of  the  torus.  This  observation  is  in  agreement  with  the  well 

known  fact  that  the  boundary  conditions  fail  to  characterize  the  solution 

17 

if  the  flow  domain  is  rultiply  connected.  It  follows  from 

'This  statement,  of  course,  in  no  way  contradicts  the  uniqueneae 
theorems  of  hydrodynamics,  which  require  also  the  specification  of  the 
initial  conditions;  see,  for  exanple,  £16^,  p.  U22. 


uniquely 


it 


4 


! 


mm 
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\ 


f 


Equation  (52)  that  the  velocity  potential  represented  by  Equation  (65)  is 
multiple -valued  unless  c is  chosen  appropriately.  In  view  of  Equation  (52), 
a circulation-free  potential  is  assured  if  and  only  if 


co 

II*- 

n = 0 


(66) 


or,  by  Equation  (63), 


-»[£  ‘JtSrfr-J 


-1 


(67) 


n = 0 


The  complete  solution  to  our  problem  is  now  given  by  Equations  (60), 

(65),  where  and  are  defined  in  Equations  (U8),  (56),  (51),  and 

the  coefficients  &n  follow  from  Equations  (63),  (6U),  (67)*  The  toroidal 

velocity  components  are  readily  obtained  from  Equations  (60),  (U8),  (12), 

with  the  aid  of  Equations  (l8)t 

00 


'or  + £ XI [2?2pn-V2  - (Un2-1)/i2pn-V2] co* not 

' n = 0 


oo 


V ^ “nPn  - 1/2  [V  8in  "*+  5 C°“  “l 

' n = 0 

Finally,  we  note  for  future  reference,  that  the  series  (60)  admits  the 


(68) 


rearrangement. 


n = 0 . m = 0 


(69) 


Solution  of  the  Torus  Problem  The  Potential-Function  Approach 

In  this  section  we  sketch  a second,  independent  approach  to  the  problem, 
in  tents  of  the  velocity  potential.  Although  this  method  of  attack  la  more 
cumbersome  than  the  one  adopted  previously,  it  should  nevertheless  prove  to 
be  instructive.  The  alternative  formulation  of  the  problem  as  a problem  of 


r*l 


J 
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Neumann,  requires  the  determination  of  a function  0(<tjp)  which  is  harmonic 
in  the  domain  0 - a - 2fT  , 0 - p - 0q}  and  obeys  the  boundary  and  regularity 
conditions, 


= 0,  0(ct,p)  -•»  0 = t as  (<t,p)  -►  (0,0) 


(70) 


where  0^  is  the  potential  of  the  uniform  velocity  field  at  infinity*  Re 


thus  put 


oo 

0(«,P)  = » ♦ ^2  bn  K 

n = I 


(71) 


with  reference  to  the  definition  (U8)  of  the  potential  functions 
Applying  the  first  of  conditions  (70),  we  reach  with  the  aid  of  Equation 
(28). 

<bn+l  - V'A*1A<V  ' <bn  “ “ 


-3/2* V “ SUlA^] 


(72) 


bQ  - 0,  (n  = Of  Ip  *••) 


Equation  (72)  is  a finite  difference  equation  of  the  second  order  for  the 
unknown  coefficients  of  superposition  bQ(n=l,  2,  3,  •••)  which  possesses 
an  integrating  factor  1/2 (**©)• 

Now  let 


cn  = ^bn"  bn  - l^Pn  - l/2Pn  - 3/2 


(73) 
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and  Equation  (72)  may  be  written  as 


argument  of  the  Legendre  functions  ip pearing  in  this  and  the 
following  equations  is  understood  to  be  qQ. 


1 1 

» I 


I 


! 


=n+  1 - '»  “ ^ [«i-  3/2  - «A*  Vi]  PA-  V2  (7U) 

This  difference  equation  of  the  first  order  Is  integrablej  by  use  of  the 
bilinear  identity  (25),  we  obtain 


22  i 


n-  1/2  n- 3/2 


= b - b , 
n n - 1 


3/2  L Zn"il  Pn-l/2 


in  which  A is  an  arbitrary  constant  of  integration*  A second  integration, 
based  on  Equation  (75),  yields 


(n  — 1,  2,  3»  •*•) 


the  notation  being  that  of  Equation  (6U)«  In  order  that  the  series  (71) 
be  convergent,  it  is  nececsary  that 

lin  bn  = 0 (77) 

n -*■  oo 

This  formal  convergence  condition,  which  here  takes  the  place  of  the  single- 
valuedness  condition  (66), ~7  serves  to  determine  the  parameter  A « Indsad, 
since 


11m  nl„  = 0 
n 

n co 


we  conclude  from  Equations  (77),  (76),  that 

A * -£*•  (79) 

2q 

with  c define * by  Equation  (67),  and  thus, 

^Recall  that  the  functions  ffn  entering  Equation  (71),  and  the  con- 
jugate stream  functions  Wn  4re  all  single-valued. 


f 


a 


1 


» 
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b 


n 


(a  = 1»  2t  3»  •••) 


(80) 


the  coefficients  eQ  retaining  their  previous  definition  given  by  Equa- 
tions (63)j  (bh),  and  (67).  This  completes  the  independent  dete  mi  nation 
of  the  velocity  potential  introduced  in  Equation  (71)* 

With  a view  toward  demonstrating  the  equivalence  of  the  solutions 
derived  in  this  and  the  preceding  sections}  we  note  from  Equation  (80)  that 


(n  — 0}  1}  2.  «••);  b^  — 0 


Consequently}  Equation  (69)  nay  be  written  as 

Pm £ * [sot  - slrt X b.]  < K-  (•» 

v. « 6 ■ = 6 

It  now  follows  free  Equations  (83)}  (U9)  that  the  conjugate  velocity  po- 
tential appears  as 


ao 


u 


n+l‘ 


(8U) 


which}  after  a permissible  rearrangement}  assumes  the  form} 

0 = a -i?  0 + "T^b  * 

T 0 Zl-*  n n 

n = 0 


(85) 


Since  = 0}  Equation  (85)  is  identical  with  Equation  (71).  It  is 
-interesting  to  observe  that  the  representation  (71)  for  in  con- 

trast to  that  given  by  Equation  (65)}  does  not  involve  any  elliptic  functions 
or  integrals. 


! 


2h 


Discussion,  Numerical  Example 

A comparison  of  Equations  (69 f or  (8h)  with  the  second  of  Equations 
($7);  or  of  Equation  (71)  with  the  first  of  Equations  (57),  reveals  that 
the  leading  term  In  the  series  contributions  to  (//' {0 Lfi)  and  0(<tjp)  are 
precisely  the  stream  and  potential  functions  belonging  to  the  ring-doublet* 
This  observation  identifies  Streeter's  solution  £ll^,  *or  ^orU8  °f  nearly 
circular  cross  section,  as  a first  approximation  to  the  exact  solution  pre- 
sented here.  The  fact  that  already  a one-tern  approximation  corresponds 
to  a torus  whose  cross  section  is  very  close  to  a circle,  is  accounted  for 
by  the  extremely  rapid  convergence  of  the  series  under  consideration* 

Although  the  convergence  of  the  representation  (69)  for  the  stream 
function  is  even  more  favorable  than  that  given  in  Equation  (60),  It  was 
still  found  to  be  more  expedient  to  use  the  latter  representation  in  the 
numerical  example  which  follows.  The  reason  for  this  lies  in  the  loss  of 
numerical  accuracy  caused  by  the  difference  formations  inherent  in  Equa- 
tion (69). 

The  stream  pattern  shown  in  Figure  2 applies  to  a shape-ratio 
A = b/a  = 3.  It  was  obtained  by  choosing  various  fixed  values  of  z,  and 
by  subsequently  selecting  <t  and  p Judiciously  and  consistent  with  the 
last  of  Equations  (7).  For  these  values  of  a and  p the  corresponding 
values  of  IfJ  were  co^mted  from  Equation  (60)  with  the  aid  of  the  numeri- 
. cal  tables  £6j.  A maximum  of  five  terms  in  the  series  for  was  needed 
in  order  to  reach  an  accuracy  of  four  significant  figures*  Once  a suf- 
ficient number  of  s-profiles  of  the  stream  surface  had  been  established, 
the  streamlines  ^-’constant  Iwrc  determined  graphically.  - -In  order 
to  locate  the  stagnation  point  — the  point  of  intersection  of  the  separa- 
tion stream ^ine  with  the  meridional  section  of  the  torus  — Equation  (68) 


25 


waa  solved  for  v^( ot>p)  = 0.  A suitable  method  of  successive  approxi- 
mations revealed  that  v^a^p^  vanishes  to  four  significant  figures  for 
cos  cl-  0.121. 
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FIG  I - TOROIDAL  COORDINATES.  TRACES  OF  SURFACES 
a = CONST  AND  £ = CONST.  ON  A HALF -PLANE 
y CONST 
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